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DIURNAL ANISOTROPY DURING SOLAR ACTIVITY CYCLE TWENTY 
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ABSTRACT 

Underground muon telescope data obtained at 
Embudo and Neutron monitor data obtained at Deep 
River are divided into two sets; one covers the 
ascending phase of the cycle (1965-70) and the 
other covers the descending phase (1971-76). The 
amplitude of diurnal anisotropy calculated from the 
data does not agree with the value predicted by the 
s imp 1 if ied version of the Diffusion-Convection 
Model (DCM) ; the discrepancy is worse for neutron 
data. 


INTRODUCTION . We report here an extension of our earlier 
study (Ahluwalia and Ericksen, 1970, 1971). Coupling 
functions (W) given by Murakami et al. (1979) and those 
given by Lockwood and Weber (1967) are used to calculate the 
amplitude of diurnal anisotropy. Median rigidity of 
response (R ) ,v 16 GV for Deep River neutrons and ~ 134 GV for 
Embudo muons. Unlike Subramanian (1971) we are not able to 
reconcile our calculations with the predictions of the 
simpl if ied version of Diffusion-Convection Model (Forman and 
Gleeson, 1975). 



DATA . Figure 1 sum- 
marizes the available 
data on the annual 

mean amplitudes of the 5 .35 
diurnal variation ob- % 
served at Embudo (open < .3oj 
circles) for the 2 
period 1965-79 and at 
Deep River (solid 

dots) for the period 

1962-79. The flags 
indicate the standard 
deviation, computed 
from the observed 

scatter in the data 
for a given year. 
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Zurich sunspot numbers _ 

(R z ) indicate the F* » G 1. 

level of solar activity. Epoch of solar polar field rever- 
sal is also shown; after 1970 the solar polar field is 
directed out of the sun in the northern hemisphere. In 
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general the amplitudes 
have the lowest values 
during activity minima 
and large values 
during activity maxi- 
mum. In particular, 
we wish to point out 
that the amplitudes at 
Embudo have large 
values for the years 
1973, 1974 and 1975 

during which the mean 
level of the solar 
activity is only about 
one third of the maxi- 
mum value in 1968. 

Bulk velocity (V) of 
solar wind is high 
during these years. 

F i G £ 

Another important feature of the dependence of the 
diurnal variation on the solar activity is shown in Figure 
It gives the annual mean time of maximum of the diurnal 
ob8erved at Embudo, Deep River, Hermanus (R a=20 

woil2w^ UanCa> T ( V 3 ° GV) ' The data e * h i b it a remarkaBle, 
worldwide, systematic, rigidity-dependent shift of the time 

Derio^^qvi 0 ^* 3101 ^ 1 v <? r i afcion to earl y hours, durinq the 
period 1971-76. The data for Huancayo and Hermanus are 

staggered m time, with respect to Deep River, to minimize 
an overlap. The rigidity dependent nature of this effect 
was emphasized by us elsewhere (Ahluwalia, 1976). Note that 
prior to 1971 the observed time of maximum at Embudo is 
nearly steady at about 15 hours (LT) but in 1976 the time of 

return' ll * a ^° ur ? (LT) * Th * diurnal times of maximum 
Aarawal values by about 1979, at all stations. 

Agrawal (1983) has reported similar behavior for neutron and 
suiface level muon detectors at sites in Canada. Also sum- 
anzed in Figure 2 is the information about the size of the 
coronal holes during solar activity cycle (SAC) 20. 

ANALYSIS we subdivide the data into two sets; one set 

covers the period 1965-70 and the other set covers the 
period 1971-76. The reader is reminded that the observed 

during 1 ^ 18 . r ather steady, on a year to year basis, 

lr ing the first period and undergoes a remarkable system- 
■atic change during the second period. The characteristics 
. _ the annual mean diurnal variation observed at Embudo and 

tL D hI? m RlVer ^ f °^ tKe tW ° P eriods are displayed in Fig. 3 on 
the harmonic dial. The error circles are for 2o. One can 

ee that the amplitudes for the two periods remain the same 
for both Embudo muons and Deep River neutrons. However, “thl 
time of maximum is significantly early for both detectors, 
for the second period (1971-76). The fact that the ampli- 
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tudes are the same for muons as 
well as neutrons, implies that 
upper cut-off rigidity R has 
the same average value for c both 
periods; in fact R ^ 70 GV. 

Also the fact that tfie vectors 
for the two periods, for both 
detectors , are well-resolved 

implies that some new effects 18 i 

contribute to the observed 
diurnal variation during 1971-76 
period . 

Diffusion-Convection Model 
(DCM) predicts that if K - 0 

and there is no net radial 
streaming (S ) and if perpendi- 
cular gradient is zero, cosmic 
rays undergo corotation. In 

this case the anisotropy ( 6 ) 12 

recorded by a detector with a 

mean asymptotic latitude X of 

viewing is given by F v G 2> 

<5 = 0.6 Cos X, % 

So we expect that 6 should be independent of solar wind bulk 
velocity (V) . An obvious limitation of the model is that it 
does not contain the concept of R . Somehow the diffusion 
coefficients Kj., K, must depend upon primary rigidity (R) 
in such a way as to lend some physical significance to the 
concept of R . The theory in its present state is obscure 
on this point. Let us assume that, 

( 0.6 cos X %, if R^R 

c 

0.0, if R>R 

c 

This form of variational spectrum is consistent with the 
results obtained in our analyses earlier, for the first 
period (1965-70). Table 1 gives a summary of results for 
the two detectors, for the three periods, assuming R =70 GV. 
The direction of anisotropy in free space is caSculated 
after applying corrections for the "orbital-effect" and for 
the "geomagnetic bending" of the primaries. These 
corrections are model independent. 

RESULTS . (1) Diurnal Anisotropy in free-space is in the 

direction of 18 hours (LT) for the first (1965-70) and third 
periods (1965-76), but the direction is significantly 
earlier for the second period (1971-76). 

(2) Within errors of observations the amplitudes of 
diurnal variations for the two detectors, for the three 
periods, are the same. The amplitude of the anisotropy, 
calculated from the data for the two detectors, is tabulated 
in the last column of Table 1. The calculated value is 
about 33% less for neutrons and about 8% less for muons when 
compared to the theoretical value. So our results do not 
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TABLE 1 




r R 

i r 

Free Space Direction 

& calc 

Station 

X 

* 

R WdR 
y o 

of Diurnal Anisotropy 


deg 

% 

_ hours (LT) 

1965-70 7¥, 


o/o 

Deep River 

00 

0.78 

18.21 ± 16.87 ± 

17.54 ± 

0.43 ± 




0 . 12 0.12 

0.12 

0.01 

Embudo 

19.1 

0.28 

17.83 ± 16.87 ± 

17.74 ± 

0.55 ± 




0.11 0.11 

0.11 

0.02 


agree with the predictions of the simplified DCM, for SAC 
20. However in practice we may not have K - 0. Also K. 
may be inversely proportional to the square of the particle 
rigidity (Jokipii, 1971). If we follow Jokipii' s suggestion 
and take (K , /K | j ) “ 0.15 and further note that median primary 
rigidity or response for muons is about a factor of eight 
larger than for neutrons, the calculated values agree much 
better with the theoretical value. But this does not 
explain why the direction of anisotropy is significantly 
earlier in the second period (1971-76). Also it is hard to 
believe that the ratio K^/K^ ^ is invariant in time. 

(3) In recent years Jokipii has emphasized very 

strongly that drifts should play an important role in the 
modulation of cosmic rays (Kota and Jokipii, 1983 and 
references therein). It would be interesting to see if one 
obtains a better agreement between theory and observations 
by invoking charged particle drifts. 
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